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Theoretical description of oscillations of electron liquid in large metallic nanospheres (with radius 
of few tens nm) is formulated within random-phase-approximation semiclassical scheme. Spectrum 
of plasmons is determined including both surface and volume type excitations. It is demonstrated 
that only surface plasmons of dipole type can be excited by homogeneous dynamical electric field. 
The Lorentz friction due to irradiation of electro-magnetic wave by plasmon oscillations is analyzed 
with respect to the sphere dimension. The resulting shift of resonance frequency turns out to 
be strongly sensitive to the sphere radius. The form of e-m response of the system of metallic 
nanospheres embedded in the dielectric medium is found. The theoretical predictions are verified 
by a measurement of extinction of light due to plasmon excitations in nanosphere colloidal water 
solutions, for Au and Ag metallic components with radius from 10 to 75 nm. Theoretical predictions 
and experiments clearly agree in the positions of surface plasmon resonances and in an emergence 
of the first volume plasmon resonance in the e-m response of the system for limiting big nanosphere 
radii, when dipole approximation is not exact. 



I. INTRODUCTION 



Experimental and theoretical investigations of plasmon excitations in metallic nanocrystals have received recently 
much attention due to possible applications in photo-voltaics and microelectronics. A significant enhancement of 
absorption of the incident light in photodiode-systems with active surface covered with metallic particles (of Au, Ag 
or Cu), with radius ten to several tens nanometers and with planar density ~ 10 8 /cm 2 , was observedi^^ 5 -^. This is 
due to mediating of light energy transfer by surface plasmon oscillations in metallic nano-components. These findings 
are of practical importance towards enhancement of solar cell efficiency especially for thin film cell technology. On 
the other side, hybridized states of surface plasmons and photons result in plasmon-polaritons^ which are of high 
significance for applications in sub-diffractional photonics and microelectronic o 9 ' 10 

For finite size crystals where the surface strongly affects the plasmon spectrum, surface plasmons occur and dominate 
the electro-magnetic (e-m) response of the metallic system. A strong dependence of resonance surface plasmon 
frequencies on the nanoparticle size and shape are reported in the literatur e) i 11 . The plasmon resonances of noble 
metals, such as gold and silver (including also copper), are of particular interest due to their frequencies located within 
the visible part of the e-m spectrum. 

Plasmon oscillations in metallic nanospheres can be excited by time-dependent electric field signal. There are 
different types of plasmon oscillations in the case of metallic nanoparticle. General types are volume and surface 
plasmons, referring to oscillations of internal electron density and surface electron density, respectively. The surface 
plasmons are linked to translational motion of all electrons which results in surface density oscillations only. The 
volume modes are related to compressional oscillations. Note that separation of these both types of collective excita- 
tions of electrons in metallic nanocrystal repeats the similar distinguishing of collective modes in atomic nucle i 12 ' 13 ! 14 
(confirmed by giant resonance experiments). For metallic ultra-small clusters with number N of electrons (in the 
range form N = 8 to N = 200 in Na clusters)^ 5 , the decoupling of volume and surface excitations is demonstrated 
by microscopic modeling^ 5 - at approximately N = 50. For ultra-small metallic clusters quantum shell effects^ 5 , and 
spill-out of electron cloud beyond the ionic jelliu m 16 i 17 disturb separate formation of surface and volume collective 
excitations, while for large clusters (with radius larger than 10 nm) the role of shells and spill-out is considerably 
reduce d 15 ' 16 ! 18 , and both modes are well defined. 

In order to excite the volume type oscillations an electric dynamical field inhomogeneous on the nanosphere scale 
is necessary while a homogeneous field excites only surface plasmons. For spherical symmetry, all modes of plasmon 
oscillations can be represented by spherical harmonics in terms of l,m, angular momentum (multiplicity) numbers. 
A dynamical electric field homogeneous over the sphere can induce only I — 1, i.e., dipole-type surface oscillations. 

The surface plasmons have been originally considered by Miei£, who provided a classical description of oscillations 
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of electrical charge on the surface of the metallic sphere within the classical model. The dipole-type Mie oscillation 
energy is not dependent of the sphere radius, in contradiction to experimental observations, both in the case of 
small and larger nanospheres. For low radius, of order of single nanometers, besides mentioned above spill-out, 
the electron-electron interactions are importan t 11 ' 1 ^ 16 ' 18 (including decay of plasmons into particle-hole pairs with 
similar energy, called as Landau dampin g 17 i 20 ) — these quantum effects influence on position of resonance frequency. 
For large spheres, with radius bigger than 10 nm, even stronger shifts of resonance are observed, probably connected 
with another mechanism, since with radius growth quantum effects turn out to be not so important as for ultra-small 
clusters. The case of bigger nanospheres is, however, of particular significance as such metallic nano-componcnts 
would be applied to enhancement of photo- voltaic effect in metallically modified solar ceUsi^^^. 

Plasma excitations in metallic clusters were analyzed within many attitudes, addressed, however, mostly to small 
clusters. In particular were developed numerical methods of calculus 'ab initio' including Kohn-Sham 'local den- 
sity approximation [LDA], similar as applied in chemistry for large molecule calculations (limited, however, to few 
hundreds of electrons) 15 i 16 i 21 , Also variational methods for energy density, semiclassical approach^ and random- 
phase-approximation (RPA) numerical summations were applied (e.g., for clusters of Na with radius ~ lnm) 22 . 
Emerging of the Mie response from the more general description was analyzed, but including only (in a numerical 
manner) single breathing volume mode^ 2 .. Commonly the 'jellium' model was applied, allowing for adiabatic approach 
to background ion system. In the 'jellium' model all the kinetics concerns electron liquid screened by uniform and 
static background of positive ion o 15 i 23 ' 24 . 

Below we present the simplified RPA-type theory of plasmon excitations in metallic nanosphere embedded in 
dielectric medium adjusted to large nanospheres (with radius above 10 nm), when quantum corrections beyond 
semiclassical approximation are not so significant as in the case of ultra-small cluster a 15 i 16 ' 18 . For optically induced 
plasma oscillations, the wavelength of incident light which excites resonance oscillations in metallic nanospheres (Au 
and Ag within several to several tens nm for radius) is considerably longer (A ~ 400 nm) than nanosphere dimension. 
Thus, the dipole-type approximation is valid, i.e., one can assume that the electric field of the incident e-m wave is 
homogeneous over the nanosphere. Therefore, only dipole type oscillations of surface plasmons will contribute the 
resonance (except for limiting big nanospheres, when also volume excitations seem to enter e-m response due to not 
exact dipole approximation). 

The experimentally observed red-shift of dipole Mie plasmon resonance for ultra-small clusters is caused mainly by 
significant spill-out effect reducing density of electron o 15 ' 16 ! 17 ' 18 . The Mie frequency is proportional to square root of 
the electron density and thus one arrives with reduced its value due to spill-out beyond the edge of ion jellium. With 
growing radius this effect weakens, as being of surface type and thus proportional to inverse radius, quite oppositely 
as red-shift observed for bigger nanospheres. For the radius over 10 nm the red-shift experimentally observed is much 
stronger than that one for small clusters and is sharply growing with radius enhancement. In order to explain this 
phenomenon observed^ in nanospheres of Au and confirmed by more precise measurements for Au and Ag nanospheres 
reported in the present paper, we have included damping of plasmon oscillations via irradiation effects, which seem 
to dominate plasmon energy losses at larger scale of radii and cause strong red-shift of the resonance. We performed 
measurements of light extinction by nanospheres (in water colloidal solution) of Au with radii from 10 nm to 75 nm 
and Ag from 10 nm to 40 nm. The resulting data reveal a strong shift of resonance towards higher wave-lengths with 
radius growth. Radiation losses, which, as we suppose, are responsible for radius dependent red-shift of resonance 
frequency, can be described in terms of Lorentz friction 2 ^, and calculated also independently in Poyting vector terms in 
a far-field zone of plasmon radiation. Damping of plasmons is caused also additionally by electron scattering processes 
and we verify that this channel is not important for a > 20 nm (a — nanosphere radius) when radiation losses are much 
stronger. The shift of the resonance frequency of dipole-type surface plasmons resulting due to damping phenomena 
is compared with the experimental data for various nanosphere radii. Emerging of the first mode of volume plasmons 
is experimentally observed for 2a ~ 150 nm for Au and 2a ~ 80 nm for Ag, due to breaking the dipole approximation 
at this range of radius, in agreement with the presented theory predictions. 

The paper is organized as follows. In the first paragraph, the RPA theor y 26 i 27 i 28 is generalized for the confined 
system of spherical shape. In the second one, the equations for volume and surface plasmons are solved (with 
particularities of calculus in the Appendix). The third paragraph contains a description of the Lorentz friction for 
surface plasmons oscillations of the dipole-type. In the fourth one, an analysis of the radiation losses is presented using 
the Poyting vector of plasmon radiation in far-field region, which supports the previous Lorentz friction estimation. 
The following paragraph comprises the summarizing of whole e-m response of metallic nanosphere system. The last 
paragraph presents a comparison of the theoretical predictions with experimental data of e-m response features of the 
colloidal water solutions with nanospheres of Au and Ag with several radii of metallic nano- components. 
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II. RPA SEMICLASSICAL APPROACH TO ELECTRON EXCITATIONS IN METALLIC 

NANOSPHERE 

A. Derivation of RPA equation for local electron density in spherical geometry 



Let us consider a metallic sphere with a radius a, located for the starting model in the vacuum, e = 1, p = 1 and 
in the presence of a dynamical electric field (magnetic field is assumed zero). The model jelliu m 15 ' 23 ^ 4 is assumed 
in order to account for the screening background of positive ions in the form of static uniformly distributed over the 
sphere positive charge: 

n e (r) = n e Q(a - r), (1) 

where n e — N e /V with n e \e\ is the averaged positive charge density, N e is the number of collective electrons in the 
sphere, V — is the sphere volume, and 8 is the Heaviside step-function. After neglecting the ion dynamics within 
the jellium model, which can be adopted in particular for description of simple metals, as noble, transition and alkali 
metals, we deal with the Hamiltonian for collective electrons, 
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n e (r )d 3 r 
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■etp{rj,t) 
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AE, 



(2) 



where rj and m are the position (with respect to the dot center) and mass of jth electron, AE represents the 
electrostatic energy contribution from the ion jellium, and ip(r, t) is the scalar potential of the external electric field. 
The corresponding electric field E(rj, t) = — Vy>(Vj, t). Assuming that space-dependence of E is weak on the scale of 
the sphere radius a (i.e the electric field is homogeneous over the sphere), (p(rj,t) = —rj ■ E(t). 
A local electron density can be written as follows^: 



p(r,t) =<V(r e ,t)\J2 S ( r - r i)\y( r e,t) >, 



(3) 



Ul dt 



He'i'^ejt), r e = (ri,r2, ...,rjf), with the Fourier picture: 



p(k,t) = J p(r,t)e- ikr d 3 r =< *'(r e ,t)\p(k)\*'(r e ,t) >, 



(4) 



where the 'operator' p(k) = ^ 



Using the above notation one can rewrite H e , in analogy to the bulk case^l, in the following form: 
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where: n e (fe) = J d 3 rn e (r)e 
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(p{k) = J d 3 rip(r,t). 
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Utilizing this form of the electron Hamiltonian one can write the motion equation for p(k): 
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or, after some algebra: 
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Jd 3 qSp(k-q)^Sp(q), 



where Sp(k) — p(k) — n e (k) is the 'operator' of local electron density fluctuations beyond the uniform distribution. 
Taking into account that: Sp(k,t) =< ^> (t)\6 p(k)\^> (t) >= p(k,t) - n e (k) we find: 
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^ / d 3 q6p(k - q, t)(k ■ q)(p{q, t) - ^ J d 3 qH < *\S P (k - q)5p(q)\* >, 
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One can simplify the above equation using the assumption that Sp(r,t) = J e lkr 5p(k,t)d 3 k only weakly varies 
on the interatomic scale, and hence three components of the first term in right-hand-side of Eq. (jHJ) can be estimated 

as: k 2 VpSp(k, i), k 3 v F /k F Sp(k, t) and k 4 v 2 7 /k 2 ^8p(k 1 1), respectively (1/fcr is Thomas-Fermi radius 2 ^, k F = !i 



e F — the Fermi energy, v F — the Fermi velocity). Thus the contribution of the second and the third components of the 
first term can be neglected in comparison to the first component. Small and thus negligible is also the last term in 
right-hand-side of Eq. {5]), as it involves a product of two 5p (which we assumed small 8p/n e << 1). This approach 
corresponds to random-phase-approximation (RPA) attitude formulated for bulk meta l 26 i 27 (note that Sp(0) — and 
the coherent RPA contribution of interaction is contained in the second term in the right-hand-side of Eq. ([8])). The 
last but one term in Eq. l[5])) can be reduced if to confine only to linear terms with respect to Sp and (p.. Next, due 



to spherical symmetry, < '5 | £ e 

j 

transform, Eq. |8]) attains the form 
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According to Thomas-Fermi approximation 2 ^ the averaged kinetic energy can be represented as follows: 
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Taking into account the above approximation and that V0(a — r) = ~^S(a 
that f(r,t) = —r ■ E(t), one can rewrite Eq. ([9]) in the following manner: 
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as well as 
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In the above formula oj p is the bulk plasmon frequency, oj 2 
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solution of Eq. (fTTj) can be decomposed into two parts with regard to the domain: 
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corresponding to the volume and surface excitations, respectively. These two parts of local electron density fluctuations 
satisfy the equations: 
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It is clear from Eq. (Ti~3"|) that the volume plasmons are independent of surface plasmons. However, surface plasmons 
can be excited by volume plasmons due to the last term in Eq. (Tj"4l (corresponding to 'a surface tail' of volume 
oscillations), which expresses a coupling between surface and volume oscillations in the metallic nanosphere within 
the above semiclassical RPA approach. 

In a dielectric medium in which the metallic sphere can be embedded, the electrons on the surface interact with 
forces e (dielectric constant) times weaker in comparison to electrons inside the sphere. To account for it, one can 
substitute Eq. (fl"4|) with the following one (Eq. (|13|) does not change): 
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Let us also assume that both volume and surface plasmon oscillations are damped with the time ratio tq which can 
be phenomenologically accounted for via the additional term, __1 dSp ^' t ' > , to the right-hand-side of above equations. 
They attain the form: 

d 2 S Pl (r,t) i 2 dS Pl (r,t) 2 e F 2 2 

and 

- l^^ 2 (r,t) + g^V/d 3 r 1F ^(^ 1 (r 1 ,t)e(a-r 1 ) (17) 
+ i5p 2 (r 1) t)6(r 1 - a)) + ^ . *( a + e _ r ). 

From Eqs (JT6J) and (fl7)) it is noticeable that the homogeneous electric field does not excite the volume-type plasmons 
but only induces the surface plasmons. 

The derived above equations for plasmon excitations in spherical metallic system are in agreement with other similar 
semiclassical approximations reviewed e.g., in Ref. [l8l 



B. Solution of RPA equations: volume and surface plasmons frequencies 

Eqs (fTBI [T7|) can be solved upon imposed the boundary and initial conditions (cf. Appendix [A"|l . Let us represent 
both parts of the electron fluctuation in the following manner: 

Spi{r,t) = n e [fi(r) + F(r,t)] , forr<a, ^ 
5p 2 {r,t) = n e f 2 (r) + a(fl,t)6(r + e- a), for r > a, (r -> a+), ^ ' 

and let us choose the convenient initial conditions, F(r, t)\t=o — 0, cr(0, i)|t=o = 0, (f2 = (9, ip) — the spherical angles), 
moreover, (l + / 1 (r))| r=0 = / 2 (r)| r=0 (continuity condition), F(r, t)\ r=a = 0, J p(r, t)d 3 r = N e (neutrality condition). 
We arrive thus with the explicit form of the solutions of Eqs (1161) and (jTTJ) (cf. Appendix |A"|) : 
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where = ^J^Ii+i/2(C) is the spherical Bessel function, Yi m (Q) is the spherical function, lu ti i = Lu p ^Jl + fc ^" a ' 2 

are the frequencies of electron volume free self-oscillations (volume plasmon frequencies) , x„i are nodes of the Bessel 

function ji (£) , loqi = uj p \J 2 i+\ are the frequencies of electron surface free self-oscillations (surface plasmon frequencies) , 

and k n i — x n i/a; lo' = . hJ 2 — \ are the shifted frequencies for all modes due to damping. The coefficients Bi m and 
V T ° 

Mmn can be determined by the initial conditions. As we have assumed that Sp(r,t = 0) = 0, we get Bi m = and 
Mmn = 0, except for I = 1 in the former case (of Bi m ), corresponding to response to homogeneous electric field. This 
mode is described by the function Qx m (t) in the general solution pTjl . The function Qx m (t) satisfies the equation: 

d 2 Qi m (t) , 2 aq lm (t) , , ,2 n (f \ 

— d&zr + n — si — h w i yitnw , 22 

f [E z (t)S m0 + V2(E x (t)S ml + E y {t)6 m -i)] , 



where uj\ = u>oi = (it is a dipole-type surface plasmon Mie frequency^). Only this function contributes the 
dynamical response to the homogeneous electric field (for the assumed initial conditions). From the above it follows 
thus that local electron density (within semiclassical RPA attitude) has the form: 



p(r,t) = p (r) +p 1 (r,t), 
with the RPA equilibrium electron distribution (correcting the uniform distribution n e ): 



Po(r) 



n e [1 + /i(r)] , for r < a, 
n e h(r), for r > a, r — > a+ 



and the nonequilibrium part, of surface plasmon oscillation type: 

0, for r < a, 



Pi{r,t) 
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(23) 



(24) 



(25) 



In general, F(r,t) (volume plasmons) and <x(f2,t) (surface plasmons) contribute to plasmon e-m response. However, 
in the case of homogeneous perturbation, only the surface I — 1 mode is excited. 

For plasmon oscillations given by Eq. ([25]) one can calculate the corresponding dipole, 

(26) 



f Air 
D(t) =e d 3 rrp(r,t) = — eq{t)a 3 , 



where Qu(t) = J^-q x {t), Qi-\{t) = \/^§-q y {t), Qw(t) — J ^-q x (t) and q(t) satisfies the equation (cf. Eq. 



-E(t). 



(27) 



III. LORENTZ FRICTION FOR NANOSPHERE PLASMONS 



The nanosphere plasmons induced by a homogeneous electric field, as described in the above paragraph, are them- 
selves a source of the e-m radiation. This radiation takes away the energy of plasmons resulting in their damping, 
which can be described as the Lorentz friction^. This damping was not included in To in Eq. ([22]) . The e-m wave 
emission which causes electron friction can be described as the additional electric field^, 



E r 



2 d 3 D{t) 



(28) 



3ev 2 dt 3 ' 

where v — is the light velocity in the dielectric medium, and D(t) is the dipole of the nanosphere. According to 
Eq. (|2"6"|) we arrive at the following relation, 



Substituting it into Eq. ([27]) we get 

dt 2 + 
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If one assumes the estimation 



d\(t) 

dt 3 



,2dq(t) 



(resulting from perturbative method of solution of the above equa- 



tion), then one can include the Lorentz friction in a renormalized damping term: 
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-E(t), 



(31) 



where (cf. Fig. 1), 
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where we used for ^ ~ + (As is the free path in bulk, vp the Fermi velocity, and C ~ 1 is a constant) 29 i 30 
which corresponds to inclusion of plasmon damping due to electron scattering on other electrons, on impurities, on 
phonons and on nanocrystal boundary. The renormalized damping causes the change in the shift of self-frequency of 

free surface plasmons, ijj' x = ^Jujf — 4^. 

Using Eq. (|32p one can determine the radius ao corresponding to a minimal damping, 

ao = — (v F c^e/2f\ (33) 



A. Radiation of surface plasmons on metallic nanosphere in far-field zone 

The Lorentz friction mechanism of energy losses of dipole surface oscillations described above can be also analyzed 
equivalently by accounting for e-m emission from oscillating dipoles of plasmons. Let us consider Eq. (|27p with 
irradiation induced damping included into r (instead of tq). For E(t) = Eq[1 — O(t)] (the rapid switching off a 
constant electric field E ) the solution of Eq. (|27[) has the form: 

/4tt en e f h for t<0, 



3 mcj^ 



e~^ T , /or i > 0, 



It is easy to calculate the loss of the total energy of the system, A = E[t = 0) — £(t = oo), i.e., by taking into 
account both kinetic and potential energy of electron system. Only potential interaction energy of oscillating electrons 
contributes, and £(t) = const. + ^a 3 q 2 (t), [the time dependent part of energy is caused by interaction of excited 
electrons, al^l J d 3 ri ,d 3 r 2 YMn^a+e^YMn^Sia+e.-r,) ^ ^ ^ ^ _> , £l > e 2 ]. For q given by Eq. flM} we 
obtain 

e 2 3 47r f en e E 



A = 8(t = 0)-S(t = oc) = -a* Y ^) . 
Radiation of the corresponding dipole, Eq. (f2"6")) . far from the sphere can be described by potentials of retarded 

r dD(t— — \ 

type, leading to the formula 2 ^ for the vector potential, A(R, t) = i — ^ t - ' . 
Hence, for far-field radiation of surface plasmon dipole oscillations we have 

B = r*A = -^nx^, (36) 



and 



E = -^=B x n, (37) 



Sill 



■e 



corresponding to the planar wave in far-field zone (n = R/ R) , with the Poyting vector II = j^E x B = 
(0 is the angle between D and R, v — cj^fe). Next, taking into account that 4£ = ^11 • ds, one can find 



lAJ oo , 2 , \ 2 

■4 = jf =3^1 ( ^ ) dt. For D 2 given by Eq. (USD, one can find the total energy transfer: 

o 6V o ^ 



en e E 
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In this way we estimated the energy loss of plasmon oscillations induced by the signal = £"(1 — 0(t)), which 
then irradiated gradually all own energy to the surrounding medium. By comparison of Eqs (135)) and l|38[) we find 




wit = 3 - — . (39) 



The above calculation of the time ratio - for oscillation damping due to radiation losses agrees with the formula 
for this parameter estimated by Lorentz friction force. 
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B. Inclusion of screening effect 



In the above consideration all irradiating electrons in the sphere were treated equivalently. Note that in surface 
plasmon oscillations take part all electrons, since it is a translational movement of all collective electrons. In fact 
some part of electron irradiation is absorbed by other electrons in the system, which reduces outside energy transfer. 
It can be accounted for in analogy to skin-effect in metals via introducing an effective radiationally active layer with 
a depth h close to the sphere surface. Thus the factor ^-(a 3 — (a — fo) 3 )/^p- can be introduced in the formula 
(|39p to account for screening skin-effect in the metallic nanosphere, with h ~ — , (a — conductivity) as for normal 
skin-cffcct 31 . Inclusion of screening results thus in reducing of Lorentz friction and in reducing of induced by radiation 
losses red-shift of resonance, from a 3 dependence to a 2 , being closer to experimental data in the latter case — cf. Figs 
2, 3 and 4. for comparison of ~ a 3 and ~ a 2 red-shift of resonance. 



IV. E-M RESPONSE OF THE SYSTEM OF METALLIC NANOSPHERES 

Let us consider number N s of identical metallic nanospheres (with radius a) randomly located in the dielectric 
medium (e > 1) of volume V. We assume the metal is simple (as considered in the previous sections) and separation 
between spheres is sufficiently large to neglect inter-sphere electric interaction. The Hamiltonian of electrons in the 

system of N s spheres (in 'jellium' model) has the form: H s = H e (r{), where ri is the position of l-th sphere 

i=i 

N s 

center, H e {ri) is the electron Hamiltonian of the Z-tli sphere. Thus the total electron wave function ip e — Yl^ e , and 

i=i 

ifc-gf- = H e {r{)w e . A density of electrons in the system has the form: p s (r, t) = ^ /°0 — r i> t): where p(r — ri,t) =< 

i=i 

^el E <K r — r i ~ r i)l^e > i s the contribution to the total electron density from the Ith sphere electrons (rj is electron 

3 

position relative to sphere center). 

The space-time Fourier picture of this electron density has the form: p s (k,u) — ^ J dtd 3 re~ lk ' r+lut p s (r ,t) = 

N s 

p(k,u>) e~ lk ri . According to the notation (|23|) we have p(k,uj) — po(k)S(uj) + p±(k,oj), with: 
i=i 



Po(k) = f d 3 rp Q (r)c 



z=i m——l 



oo Z oo 



y ' l=lm=-ln=l I p 

+n e J dnr\ J ^e-^Mkmr^YinA *("+"n,)-i("-"».) ) 



o ) 

here a = afi, r = ^. If now one uses the continuity equation, = divj s , or in the Fourier form, fc • j s (k,u>) = 

uip s (k,ui) (here j s is the electron current), then one can find: k ■ j s (k,u) = ojp\{k,oS) e~ zk ' Vl . 

i=i 

For long wave-length limit (ka << 1, which is appropriate for the e-m response of nanospheres; i.e. assuming the 
'dipolc approximation', when only linear in k terms remain) we use the following approximations: e~ r ~ 1 — ik ■ r 
and k ■ j(fc, to) ~ k ■ j(0, u>). Thus one can rewrite the continuity equation in the form: 



fe- j a (0,u) = N S U J £ E E ^r, 

Z=l m=—l n=l 



jgfr J dnrfjiikniri) J dfiiY; m (fii)(— ik ■ n) 



^ ^-S+i)^ J ^i^V^K^) / ^^(fix) (-ifc • a)] ^+^)-^-^) (41) 

OO I . , I 



£=1 m=—l 

N. 



(as for sufficiently dense nanocomponent system e n — 1j an d /<inij m (f2) = 0, for Z > 

i=i 
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1). Assuming a rapid excitation of all frequencies (E(t) = Eo2ir5(t), and thus E(ui) = Eg) 
one can write from the Ohm law — ej s (0,u>) = (t(lu)Eo, and next, for Eq in z-th direction, 
k ■ E a = cos®, fe n = fcri^|coseFio(Oi) + ^=sm(e)[Fn(ni) + yi_i(ni)]|, and k ■ a = 

fca^|cos6Fi (O / 1 ) + ^ I sm(e)[Fii(n / 1 ) + Fi_i(r2 / 1 )]|, where ^ = (6>i,0i) and Q[ = (6>i,</>i). Thus one can ob- 



tain: 



f oo a 



a(uj)E kcosQ = J if 1 E n e / i"idr 1 j 1 (k ln r 1 ) 

v U=i o 

(i + ^fc) ([^io«»0 + sin6^l(A n + Ai_i)](6(u - w nl ) - <5(w + w nl ))) ( 42 ) 
+a[Bi cos6> + sin<d^(Bn + - w i) - <5(w + w i))| , 

from which An = 5n+Bi_i = 0, while Ai = ^-Eo and B w = vn e a 3 E (the constants /i, v will be determined 

later). From the above it follows: 

*(u>) = l -^Jfn e a 4 {v(S(cj - w i) - <5(^ + woi)) 

oo x„i ^ (43) 

+ M E (24^+33.^)3.^ / x 3 dxji(x)((5(w - w„i) - S(uj + uj nl )) I , 

where u; 2 ^ = oj 2 {1 + x 2 nl / (kxa 2 )). Via the formula for one can now derive the dielectric response function e(w) = 



+00 



£'(w) + *e"(w), with e"(w) = ^cr(w) and e'(w) = £ + ±P / (to- 

— CO 



3 



+A 4 E (2fc T a2+L^ SF / dxX 3 j!(x)(S(uj-L; nl )-S(LJ+UJ nl )) 
n=l nl nl 



(44) 



and 



e (w) -£ + — ^ — V"3" r <-^ 2 

00 Xnl *) (45) 

+M E ^aW^, 1 dxx 2 n (x)V^^ . 

n=l 1,1 1,1 nl J 

00 2 2 

One can determine now the constants fi, v using the sum rule: J du)u)e"(u>) = " 27r e , where n = ^ = ^y "' , (here, 



V — the volume of the whole system, Vq — the volume of the single nanosphcre), and the static value of the dielectric 

response of the system e(0) = e + 8 J duoV-^r- (assumed to be known). These conditions give: 



<4<* / £(0)-£\ u 2 p c / £ (0)- e \ 
A* = r T 1 5 ' v = r T aa ' ( 46 ) 

00 x n \ oo a^Til 

where a x = E (2fc^WI)^ / dxx 3 ji (x) ^ , a 2 = £ (2fc . aa+ 3 3 ^ )g . / <te 3 ji(a;)^, 

n=l T nl 1,1 n=l T 711 nl 

r - N Y<1 n, - 4|eKjV 3 a 4 [77 
t-0 — Jv s v i « — ( 27r )3 y 3 v 

By introducing the oscillator strength f(u>) = 2l ^ e Ca ^ , one can express the dipole-type dielectric response of the 
considered metallically nanomodified system in a more conventional form: 

1 



£H=£V) +i £>)=£ + C0 ^]r ^ 



n 2W„1 
n— 



+ 



w„i — u) — ie u) n i + u + le 
where e = 0+, f(uj) = E /n[*(w - u nl ) + S(w + u nl )], f = (£g=£° ai - a 2 ), /„ 

^ (i - £j ^) m^h^r^ S dxx^Mx), j = e u = i. 



(47) 
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We can include attenuation (also due to irradiation losses), of dipole excitations (I — 1), via the damping term 
7 m oscillator-type equation for plasmons (which can be added to left-hand-side of Eq. (fT5)) and assuming that 

all modes are damped with the same attenuation time r) . Thus the time-dependent solution of such modified equation 

attains the form as given by Eq. (|20[) with the factor e~*/ r , and with the shifted frequency u>' n — \ lo\ \. Similarly 



to the equation for the surface plasmons, Eq. (TTi|) , can be added (to its left-hand-side) the damping term - 



2 dp 2 {r,t) 



at 



. It 



leads to the factor e */ T for the first part of the solution (f2Tj) (and simultaneously shifted frequency uj' iq = y uif Q — 

and the second term of Eq. (|2"Tj) acquires the additional factor e~ t l T (and shifted frequency u>' n — y*u^" — X, ) . The 
corresponding change in the dipole- type (1 = 1) e-m response function (|47|) resolves thus to the following expression: 



e(uo) = e'(lj) + ie"(uo) — e + c Q u, 
The above equation can be rewritten as follows: 

fn 



In 



n=0 nl 



■ c uj 



2sr 



P Z^i 0, ,' 



i - U) - - 



■ iO 



(48) 



(49) 



and 



e"{u) = c id 



fn 



(50) 



V. MEASUREMENT OF THE DIPOLE SURFACE PLASMON FREQUENCIES IN NANOPARTICLES 

WITH VARIATION OF THEIR RADIUS 

To determine the plasmon frequencies in metal nanoparticles as a function of the nanoparticle radius, extinction 
spectra of colloidal solution of Au and Ag nanoparticles with radii ranging form 10 nm to 75 nm for Au, and from 
10 nm to 40 nm for Ag, respectively, were measured. The nanoparticles, prepared as a water colloidal solution with 
an average size distribution not exceeding 8% and an almost constant total mass per ml independent of the particle 
radii, were obtained from British Biocell International. The particular data of the Au nanoparticles are listed in Tab. 
1. 



Tab. 1. Nanoparticle data for Au colloidal solutions 



nominal nanosphere radius 
[nm] 


10 


15. 


20 


25 


30 


40 


50 


75 


average nanosphere radius 
[nm] 


10.2 


15.55 


20.55 


24.65 


29.35 


39 


49.45 


77.15 


particle density [10 9 /ml] 


700 


200 


90 


45 


26 


11 


5.6 


1.7 


total volume [10 15 nm 3 /ml] 


3.11 


3.15 


3.27 


2.82 


2.75 


2.73 


2.84 


3.27 



The extinction spectra were measured using a xenon lamp operated at 150 W in combination with a Monospek 
1000 monochrometer, providing monochromatic light at wavelengths from 300 nm to 900 nm. To accommodate 
for the response of the aqueous environment, the cuvette containing the colloidal solution and the lamp spectrum, 
reference measurements were performed to which the extinction spectra are normalized. The extinction coefficient 
shown in Figs. [2] and [3] is defined as the fraction Tc /Ta u ,Acj, where Tq is the light intensity transmitted through the 
cuvette containing de-ionized water and TA u ,Ag light intensity transmitted through the cuvette containing Au or Ag 
nanoparticle colloidal solutions. 

The results are presented in Fig. [2]for Au and in Fig. [3]for Ag, respectively. The red-shift of the resonant frequency 
with growth of nanosphere radius is clearly noticeable. This is accompanied by the broadening of the attenuation 
peak and variation of peak height (at the beginning growth and next lowering of the peak height). These features are 
collected in the Fig. 2 (bottom) for Au, where the position of the center of extinction peak, its half-width and height 
are plotted versus the nanosphere radius. 

For nanoparticles of gold, silver and copper in the air, in water and in a colloidal solution, one can find ao ~ 10 — 14 
nm (cf. Eq. (|33[> . Fig. 1 and Fig. 4), i.e., the radius of nanosphere corresponding to minimal damping, which well 
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corresponds to experimental dat a 32 ' 33 . It is a cross-over point for the resonance red-shift versus a. For a > ao damping 
increases due to Lorentz friction (proportionally to a 3 , or after inclusion of screening, proportionally to a 2 ) but for 
a < ao damping due to electron scattering dominates and causes opposite behavior — enhancement of damping with 
lowering radius (proportional to -, in agreement with experimental observation s 29 ' 32 ), which leads to cross-over of 
resonance red-shift dependence on a. 

Surface plasmon oscillations cause attenuation of the incident e-m radiation where the maximum of attenuation is 

at the resonant frequency 28 uj\ — \ju>\ — This frequency diminishes with growth of a, for a > ao according to Eq. 

(|32[) . which agrees well with the experimental measurements for Au and Ag presented in Figs 2 and 3 and in Tabs 2 
and 3 (after inclusion of screening via skin-effect type correction — cf. Fig. 4). 

Tab. 2. Resonant frequency for e-m wave attenuation in Au nanospheres 



radius of nanospheres [nm] 


10 


15 


20 


25 


30 


40 


50 


75 


hio\ (experiment) [eV] 


2.371 


2.362 


2.357 


2.340 


2.316 


2.248 


2.172 


1.895 


huj^ (theory, h — a) [eV], no = 1.4 


3.72 


3.716 


3.71 


3.69 


3.66 


2.41 


2.37 


XXX 


Tiu)^ (theory, h = a) [eV] , n = 2 


2.601 


2.60 


2.59 


2.58 


2.56 


2.38 


1.64 


XXX 


Tiuji (theory, h = 6 nm) [eV], no = 2 


2.601 


2.600 


2.599 


2.595 


2.58 


2.55 


2.47 


1.84 



(C = 2, v F = 1.396 • 10 6 m/s, X B = 5.3 • 10~ 8 m, [cf. Eq. ((32 
XXX — overdamped oscillations; 



],u p = 1.371 • 10 16 1/s, wi = 3.96 • 10 15 1/s [for n=2]) 



Tab. 3. Resonant frequency for e-m wave attenuation in Ag nanospheres 



radius of nanospheres [nm] 


10 


20 


30 


40 


Tvjj't (experiment) [eV] 


3.024 


2.911 


2.633 


2.385 


Tiloi (theory, h — a) [eV], no = 1.4 


3.71 


2.70 


2.66 


2.41 


fud'i (theory, h — a) [eV], no = 2 


2.61 


2.60 


2.56 


2.382 


fiLdi (theory, h = 8 nm) [eV], no = 2 


2.60 


2.59 


2.58 


2.55 


(C=2, v F = 1-393 • 10 6 m/s, X B = 5.3 • 10~ 8 m, [cf. Eq. 132]) 


, lu p = 1.37- 10 16 



The observed behavior well corresponds with the formula 



1/s, uji = 7.89 • 10 15 1/s [for n=l]) 
with the additional skin-effect factor in the last term 



-j 1 (a 3 — (a — h) 3 )/^-a 3 (reducing a 3 to a 2 radius dependence), which gives the damping rate - for surface dipole 

l/r 2 

plasmons versus a. This damping leads to Lorentzian shape of attenuation peak (in response function), ^_ u /yi + i/ T a 

with central position at frequency u>' = y — 1/t 2 . As l/r scales as a 2 , after inclusion of screening — cf. Fig. 4, it is 
dominating contribution to the value given by Eq. (|32p for a > 20 nm (in agreement with the experiment). Inclusion 
of screening (Fig. 4) reduces irradiation-induced damping of plasmons at limiting large values of nanosphere radius 
and allows to avoid overdamped regime in this case, which is entered by unscreened a 3 damping rate dependence 
at a ~ 75 nm (cf. Tab. 2). The agreement between the model and the experimental data suggests that for the a 
dependent red-shift of resonance frequency of surface plasmons in metallic nanoparticles of large size (10 < a < 75 
nm) responsible are plasmon energy losses caused by Lorentz friction. 



VI. CONCLUSIONS 



We have analyzed red-shift of Mie frequency of dipole surface plasmon oscillations in large metallic nanospheres, 
with radius beyond 10 nm. For this region of metallic cluster size the dominating channel of plasmon damping starts 
to be radiation loss due to Lorentz friction. At approximately 10 nm for nanosphere radius the cross-over point 
of red-shift versus nanosphere radius occurs. For lower radii the j l rule dominates describing scattering Fermi type 
mechanisms of damping, for higher radii the damping due to Lorentz friction, with radius dependence ~ a 3 (or ~ a 2 
when screening is included), prevails and quickly completely dominates plasmon attenuation. The resulting red- 
shift of damped harmonic oscillation well reproduces the experimentally observed resonance positions with respect 
to metallic sphere size. We have measured the resonance positions via observation of extinction of light in water 
colloidal solutions of Au nanospheres with radii between 10 and 75 nm, and Ag with radii between 10 and 40 nm. The 
theoretical predictions well fit to the experimental behavior, especially if include corrections due to radiation screening 
of skin-effect type. The resulting a 2 scaling well reproduces the experimental curves for skin-depth of order of 6 nm 
(for Au). At the limiting value of metallic nanosphere radius (a > 75 nm) an almost overdamped oscillation regime 
is expected from theoretical analysis, earlier for Ag than for Au, due to bigger conductivity in Au and thus stronger 
reducing damping than in Ag. In experiment the corresponding large red-shift in the almost overdamped regime is 
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observed for both Aii and Ag. Moreover, the first volume mode with original energy above fao v (bulk volume plasmon 
frequency) in the case of an almost overdamped regime is strongly red-shifted and emerges in extinction features as 
an additional smaller peak on the left side of surface plasmon peak for sufficiently large nanospheres, when dipole 
approximation is not exact. 
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APPENDIX A: ANALYTICAL SOLUTION OF PLASMON EQUATIONS FOR THE NANOSPHERE 

Let us solve first the Eq. (|16p , assuming the solution in the form: 

6p 1 (r,t)=n e [f 1 (r)+F(r,t)], forr<a. (Al) 

Eq. P^|) resolves thus into: 



V 2 /iM-^/i(r) -0, 
The solution for function fi(r) (nonsingular at r — 0) has thus the form: 



(A2) 



/l(r) = a£ fc^T ( e ~ feTr " ^ ' (A3) 



where a is a constant, kx — \J 67T "^ e2 = {kr is the inverse Thomas-Fermi radius), uo p = %J 47r ^ e2 (bulk 

plasmon frequency). 

Since we assumed F(r,0) = 0, then for function F(r,t) the solution can be taken as, 

F(r, t) = F u (r)sin{uj't)e- Tot (A4) 



where lo' = \Joj 2 + 1/t, 2 . F u (r) satisfies the equation (Helmholtz equation): 

V 2 FL(r) + fc 2 F w (r) = 0, (A5) 



with k 2 = - 5 !l p . A solution of the above equation, nonsingular at r — 0, is as follows: 

F u (r) = Aj t (kr)Yi m (n), (A6) 

where A is a constant, — \A7(2£)-^+i/2(£) the spherical Bessel function [I n (Q the Bessel function of the first 
order], and Yi m {Vt) the spherical function (O the spherical angle). Owing to the semiclassical boundary condition, 
F(r,t)\ r=a — 0, one has to demand ji(ka) = 0, which leads to the discrete values of k = k n i = x n i/a, (where 
x n i, n = 1, 2, 3..., are nodes of ji), and next to the discretization of self-frequencies: 

^="K i+ |?)- (A7) 

The general solution for F(r,t) attains thus the form 

oo I oo 

F(r,t)=J2 ^^™J/(fcni0^m(^)^(^t)e- rot . (A8) 

1=0 m=—l n=l 

A solution of Eq. (TT7]) we represent as: 

Sp 2 (r, t) = n e f 2 (r) + a(Q, t)S(r + e - a), for r > a, (r -> a+, i.e.e -> 0). (A9) 
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The neutrality condition, J p(r, t)d 3 r — N e , with Sp2(r, t) = <r(u>, t)S(a + e — r) +n e /2(r), ( e — > 0), can be rewritten 

a oo a 

as follows: — J drr 2 fi(r) = J drr 2 f2{r), J d 3 rF(r,t) = 0, J dO,a(il,t) — 0. Taking into account also the continuity 
o a o 

condition on the surface, 1 + fi(a) = /2(a), one can obtain: /z(r) = f3e~ kT ( r ~ a * > / (kxr) and it is possible to fit a (cf. 
Eq. (|A3])) and (3 constants: a = ^§±±, (3 = k T a- ^§±± (l - e~ 2kTa ), which gives Eqs d 



From the condition J d 3 rF(r,t) = and from Eq. (|A8j) it follows that ^4oo« = 0, (because of / d£lYi m (oj) = 


47r<5/ ^mo)- 

In order to remove the Dirac delta functions we integrate both sides of the Eq. (JT7J) with respect to the radius 

00 

length (J r 2 dr...) and then we take the limit to the sphere surface, e — > 0. It results in the following equation for 


surface plasmons: 

I=Om=-l 

00 1 00 a r l+2 (A10) 

+Wpn e X! E A imn^^Yi m (n) J drx-^+jj^knir^siniujnit), 

[=0m=-[n=l 

+ff [^.(*)iio(n) + V2f? x (t)yu(n) + ^^y^n)] , 

where = c^p^ITi' ^ n derivation of the above equation the following formulae were exploited, (for a < r\): 
d 1 d °° a 1 °° la 1 ' 1 

75 — 7=5 20 = = 7T E — P '( cos 7) = 51 ^+T p '(cos7), (All) 

1 

where Pi(cosj) is the Legendre polynomial [P^cosj) = ^jfrr E YlmityY^^i)], 7 is an angle between vectors 

m— — I 

a = ar and r%, and (for a > ri): 

d 1 d °° r l 00 ' l + l r l 

JT , l 2 ^= = 7T E ^T^( cos 7) = - E E 4 ^^7-T^ yi ™( r! )^( fi i)- ( Al2 ) 

oa ^/ a 2 + r 2 _ 2aricos7 <?a ^ a i+i 2/ + 1 a'+^ 

Taking into account the spherical symmetry, one can assume the solution of the Eq. (|A10|) in the form: 

00 1 

<r(Q,f)=E E 1i m (t)Y lm (n). (A13) 

1=0 m=-l 

From the condition J a(cot)dD, — it follows that goo = 0- Taking into account the initial condition cr{u), 0) = we 
get (for I > 1), 

qimit) = ^fsmf^^e-'^tl - S n ) + Qi m {t)5 n 

+ E A imn to 2 -i2i > +r )w i n e J dn ^trsji (k n ir 1 )sin(u' nl t)e- t / To , ^ A14 ^ 



where u>' Ql = y — 1 Jtq and Qi m (t) satisfies the equation: 

+ ^^g 1 ^ + «&Qim(«) = ^ [E z (t)6 m0 + V2E x (t)S ml + V2E y (t)S m ^] . (A15) 
Thus cr(w,t) attains the form: 

00 z 1 
o-(fi,i) = £ E lim^^fsml^e-^H E QimW^im(^) 

l=2m=-l m=-l (A16) 

00 Z 00 (Z+l)w 2 r ! + 2 

+ E E E A nim lu2 , 2l+ ^i ym(^)n e /dri^j;(fc n ;ri)sm(^,;^)e~ t/ro - 
i=lm=-|tj=l " v ; »« 
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FIG. 1: Damping rate •=• of surface plasmons calculated according to the formula (|32[) versus nanosphere radius for Au (in 
water colloidal solution with refraction factor n = 2 in Eq. J32} C = 2, As = 5.3 • 10" 8 m), dashed line— 10 13 1/s level 



Extinction spectra of Au nanoparticles 
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FIG. 2: The results of measurement of light extinction in water colloidal solution of Au nanoparticles with radii a indicated 
in the inset; bottom — extracted red-shift of resonance frequency (left), half- width of the attenuation peak (central) and peak 
height (right) versus nanosphere radius a 
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Extinction spectra of Ag nanoparticles 
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FIG. 3: The results of measurement of light extinction in water colloidal solution of Ag nanoparticles with radii a indicated in 
the inset 
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FIG. 4: Damping rate (Eq. ([32])) of surface plasmons and corresponding resonance shift with respect to nanosphere radius for 
Au in colloidal water solution with inclusion of screening via skin layer with the depth h = 6 nm (upper) and without screening 
(lower); screening (skin-effect) improves fitting with the experimental data — cf. Fig. 2 



